Switching of discrete optical solitons in engineered waveguide arrays 
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We demonstrate simple methods for controlling nonlinear switching of discrete solitons in arrays 
of weakly coupled optical waveguides, for both cubic and quadratic nonlinearities. Based on the 
effective discrete nonlinear equations describing the waveguide arrays in the tight-binding approx- 
imation, we develop the concept of the array engineering by means of a step-like variation of the 
waveguide coupling. We demonstrate the digitized switching of a narrow input beam for up to 
eleven neighboring waveguides, in the case of the cubic nonlinearity, and up to ten waveguides, 
in the case of quadratic nonlinearity. We discuss our predictions in terms of the physics of the 
engineered Peierls-Nabarro (PN) potential experienced by strongly localized nonlinear modes in a 
lattice, and calculate the PN potential for the quadratic nonlinear array for the first time. We also 
confirm our concept and major findings for a full-scaled continuous model and realistic parameters, 
by means of the beam propagation method. 
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I. INTRODUCTION 

Discrete nonlinear systems are known to support self- 
localized modes that exist due to an interplay between 
a coupling between the lattice sites and nonlinearity 0. 
Such spatially localized modes of discrete nonlinear lat- 
tices existing without defects are known as discrete soli- 
tons or intrinsic localized modes; they appear in many 
diverse areas of physics such as biophysics, nonlinear op- 
tics, and sohd state physics 0,I3- More recently, such 
modes have been predicted in the studies of the Bose- 
Einstein condensates in optical lattices and photonic- 
crystal waveguides and circuits Q . 

One of the most important applications of discrete 
solitons is found in nonlinear optics where discrete 
optical solitons were first su^ggested theoretically by 
Christodoulides and Joseph 6] for an array of weakly 
coupled optical waveguides. Because the use of discrete 
solitons promises an efficient way to realize and control 
multi-port nonlinear switching in systems of many cou- 
pled waveguides, this field has been explored extensively 
during last ten years in a number of theoretical papers 
(see, e.g., Refs. 0, H, I3i as an example). More im- 
portantly, the discrete solitons have also been generated 
experimentally in fabricated periodic waveguide struc- 
tures (see, e.g., some original papers reporting on the 
experimental observations 0, lTl| and also the recent 
review papers 0,0,01). 

The majority of theoretical studies conducted so far is 
devoted to the analysis of different types of stationary 
localized modes in discrete nonlinear models and their 
stability. Consequently, experimental papers have re- 
ported on the observation of self-trapped states in the 
periodic systems with broken translational symmetry and 
some of their properties, in both focusing and defocus- 
mg regimes However, only very few studies 

and experimental demonstrations addressed more specific 



properties of localized modes introduced by discreteness 
such as the soliton steering in and discreteness-induced 
trapping by the lattice (see, e.g., Ref. 0). As a result, 
a very little effort has been made so far to link these 
findings with realistic applications of discrete solitons for 
multi-port all-optical switching. 

Indeed, one of the major problems for achieving con- 
trollable multi-port all-optical switching of discrete soli- 
tons in waveguide arrays is the existence of an effective 
periodic Peierls-Nabarro (PN) potential which appears 
due to the lattice discreteness. As a consequence of this 
potential, a narrow large-amplitude discrete soliton does 
not propagate freely in the lattice and, instead, it be- 
comes trapped by the array. Several ideas to exploit 
the discreteness properties of the arra y fo r digitized all- 
optical switching have been suggested Ha, 113 ■ However, 
the main result of those earlier studies is the observa- 
tion that the discrete solitons can be well controlled only 
in the limit of broad beams whereas the soliton dynam- 
ics in highly discrete arrays has been shown to be more 
complicated and even chaotic 0|. 

In this paper, we explore in detail an effective way to 
control nonlinear switching of discrete solitons in arrays 
of weakly coupled optical waveguides earlier suggested in 
our brief letter [l^- First, using the discrete model valid 
in the tight-binding approximation, we estimate the PN 
potential experienced by a strongly localized nonlinear 
mode that is kicked initially in a cubic nonlinear waveg- 
uide array. The result suggests a possible control mech- 
anism for the switching of strongly localized excitations 
by means of a step-like variation of the waveguide cou- 
pling. For particular types of the engineered arrays, we 
are able to demonstrate the digitized switching of a nar- 
row input beam for up to eleven waveguides. Second, 
we demonstrate the validity of the predictions made in 
the framework of the discrete model by performing a full- 
scaled continuous simulation using realistic parameters. 
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Last but not least, we extend the concept of control- 
lable digitized switching of discrete optical solitons to the 
case of quadratic nonlinear waveguide arrays, where the 
experimental observation of discrete optical solitons has 
been reported very recently p^ . Here, we obtain, for the 
first time to our knowledge, the PN potential for the dis- 
crete soliton and demonstrate numerically the digitized 
switching for up to ten waveguides. 

The paper is organized as following. In Sec. II we 
study the arrays of cubic nonlinear waveguides. First, 
we consider the system dynamics described by the dis- 
crete nonlinear Schrodinger equation, and show how to 
modulate the waveguide coupling in order to suppress 
the chaotic dynamics and achieve fully controllable dig- 
itized switching. We also employ the beam propagation 
method and simulate numerically a more realistic contin- 
uous model of the waveguide arrays with realistic param- 
eters, and confirm that our concept can be very useful for 
optimization of the soliton switching in realistic settings. 
Next, in Sec. Ill we extend our analysis to the arrays of 
weakly coupled quadratic nonlinear waveguides, where 
discrete quadratic solitons are composed of the coupled 
beams of the fundamental and second-harmonic fields. 
Finally, Sec. IV concludes the paper. 

II. CUBIC NONLINEAR WAVEGUIDES 

The most common theoretical approach to study the 
discrete optical solitons in arrays of weakly coupled opti- 
cal waveguides is based on the decomposition of the elec- 
tric field of the periodic photonic structure into a sum 
of weakly coupled fundamental modes excited in each 
waveguide of the array; in solid-state physics this ap- 
proach is known as the tight-binding approximation. Ac- 
cording to this approach, the wave dynamics is described 
by an effective discrete nonlinear Schrodinger (DNLS) 
equation that possesses spatially localized stationary so- 
lutions in the form of discrete localized modes. Many 
properties of the discrete optical solitons can be analyzed 
in the framework of this approach and the DNLS equa- 
tion nil. 

A. Discrete model 

1. Homogeneous Arrays 

A standard model of a weakly coupled array of cubic 
nonlinear waveguides is described by the DNLS equa- 
tion that we write in the normalized form pol | , 

^~dz ^ + ""-l) + 7l"n| Mn = 0, (1) 

where u„ is the effective envelope of the electric field in 
the n-th waveguide, the normalized parameter Y is pro- 
portional to the propagation constant of a single waveg- 
uide and it characterizes the coupling between the neigh- 
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FIG. 1: Example of a homogeneous waveguide array and the 
generation of a discrete soliton by exciting a single waveguide. 

boring waveguides, and z is the propagation distance 
along the waveguide. The parameter 7 = uj^n^j (cA^fi) is 
the effective waveguide nonlinearity associated with the 
Kerr nonlinearity of the core material. Figure^ shows a 
typical experimental structure of a quasi-one-dimensional 
homogeneous waveguide array and the excitation scheme 
for generating a discrete optical soliton. 

Steering and trapping of discrete optical solitons have 
been analyzed in the framework of the model in a 
number of theoretical studies. Being kicked by an ex- 
ternal force, the discrete soliton propagates through the 
lattice for some distance, but then it gets trapped by 
the lattice due to the discreteness effects. For a stronger 
kick, the output soliton position fluctuates between two 
(or more) neighboring waveguides making the switching 
uncontrollable 17]. 

In order to show this feature, first we consider homo- 
geneous arrays and select an input profile in the form of 
a narrow sech-like beam localized on a few waveguides, 

u„(0) = A sech[A(n - n,)/\/2] g-^'^C"-"-), (2) 

for n — ric — 0,±1, and m„(0) — 0, otherwise. For the 
particular results presented below, we consider an array 
of 101 waveguides and place the beam at the middle po- 
sition, ric = 50. The maximum normalized propagation 
distance used in our simulations is Zmax = 45. Parame- 
ter k in the ansatz Q has the meaning of the transverse 
steering velocity of the beam, in analogy with the contin- 
uous approximation. It describes the value of an effective 
kick of the beam in the transversal direction at the input, 
in order to achieve the beam motion and shift into one 
of the neighboring (or other desired) waveguide outputs. 

In our simulations, we control the numerical accuracy 
by monitoring the two conserved quantities of model fQ), 
the soliton power 

p = 5]K(^)p, (3) 

n 

and the system Hamiltonian, 

H = - E {^("«<+i + ">"+i) + (7/2) k« I'} . (4) 

n 
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FIG. 2; An example of erratic switching of a localized in- 
put beam with a slight variation of the beam intensity in a 
homogeneous array. 



The input condition does not correspond to an ex- 
act stationary solution of the discrete equation |^ even 
for fc = and, as the input kick (fc ^ 0) forces the lo- 
calized wave move to the right (fc < 0) or left (fc > 0), 
its motion is accompanied by some radiation. The effec- 
tive lattice discreteness can be attributed to an effective 
periodic potential, the PN potential, which is dynamic 
and changes in time. Due to both the strong radiation 
and the presence of the PN barrier which should be over- 
taken in order to move the beam transversally, the dis- 
crete soliton gets trapped at one of the waveguides in the 
array. In most of the cases, the shift of the beam posi- 
tion to the neighboring waveguide is easy to achieve, as 
shown in many studies jl7l |. However, the soliton switch- 
ing becomes rather complicated and even chaotic. This 
is shown in Fig. [21 where, for a fixed value of the input 
angle, a slight variation in the beam intensity results in 
a erratic switching of the beam. 



2. Modulated Arrays 

In this paper, we suggest to modulate the coupling in 
the waveguide array in order to achieve a controllable 
output and to engineer the switching results. What this 
modulation of the couplings does is to affect the PN bar- 
rier, providing us with a simple physical mechanism for 
fine tuning and control of the beam self-trapping. 

To justify the validity of our concept, we perform 
a qualitative estimate of the PN barrier in the frame- 
work of the applicability of the discrete model and per- 
turbation theory. We study the case of strongly local- 
ized modes |2lj propagating in a homogeneous waveguide 
array with identical coupling between the neighboring 
waveguides, described by Eq. QJ. We consider a gen- 
eral localized mode that we want to propagate through- 
out the array. Due to discreteness, our system lacks the 
translational invariance and, as a result, some of its en- 
ergy must be supplied in order to force the mode moving. 
Another way to look at this problem is to consider that, 
because of the lattice discreteness, the localized mode 



"sees" a potential barrier (the PN barrier), whose height 
depends on the effective discreteness of the system as 
seen by the excitation Thus, for wide modes, the 

barrier will be smaller that for narrow modes. A rough 
estimate of this PN barrier can be obtained by equating 
it to the difference in the values of the Hamiltonian, be- 
tween the mode centered at a waveguide (odd mode) and 
the mode centered between two neighboring waveguides 
(even mode) . 

In order to evaluate a change of the PN barrier for the 
mode initially kicked by an external force, we introduce 
an initial phase tilt that is proportional to the factor 
~ exp(— zfcn) in the discrete case. Our purpose is not only 
provide an extension to the earlier results [23|, but also 
study, for the first time to our knowledge, the variation of 
the effective PN potential for an initially kicked localized 
mode. 

Odd modes. We consider a strongly localized mode 
(SLM) in the form of three excited cites. 



U„{z) ^ Un e" 

« uo {0..,0,eie*'=,l,eie~ 



(5) 



where uq is the mode amplitude, fc is the parameter of the 
initial "kick" (an effective transverse angle) applied to the 
mode, Ai is the longitudinal propagation constant, and 
ei is a small parameter, to be determined from Eq. 
After substituting Eq. lO into Eq. |^ and keeping only 
linear terms in ei, we obtain 

Ai = 2eiy cos(fc) -I- 7Uq 

and ei = ycos(fc)/Ai, so that 

Ai « 7^0, ei « y cos(fc)/7UQ < 1. (6) 

Even modes. In this case, the SLM mode has the form 



« Uo {0..,0,e2e*^l,e' 



,£26 



0,...0} 



(7) 



where, as above, uq is the amplitude of the even mode, 
fc is the initial angle or effective parameter of the initial 
"kick" , A2 is the longitudinal propagation constant of the 
even mode, and £2 is a small parameter. After substitut- 
ing Eq. Q into Eq. and keeping only linear terms in 
£2, we obtain 

A2 = (1 + 62) cos(fc) 7M0 
and £2 — V cos{k)/ X2, so that 

y cos(fc) 



A2 



V^cos(fc) -I- 7U0, £2 



y cos(fc) -I- 7-Uq 



« 1. (8) 



^From Eqs. © and ^ we come to the conclusion that, 
in order to have strongly localized modes, the nonlinear 
contribution described by the term 7U§( or juq) must 
be much larger than the linear term described by the 
term Fcos(fc). Now, for calculating the PN barrier, we 
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should relate the amplitudes of the modes of two different 
symmetries. One way is to think of both the modes as 
different states of a single effective mode shifted by a half 
lattice site along the chain. This means that the power 
content of both modes must be identical, since the power 
P — J2n is ^ conserved quantity. To the first 

order in ei and £2, we obtain 

Podd = ul + O{el), Pevcn - + 0(e2) (9) 

Thus, the relation Podd = Pcvcn, implies Uq ~ 2u§. We 
are now in position to compute i?odd and i?cvon for a 
strongly localized mode, using the above relation and 
Eqs. O, ®, and ®, 

HoM « -| 4 + 0{el) 

i?cven ~ -^ut-2ulVcos{k)+0{ei-e2), (10) 

which implies that the PN barrier A'^'^^ for the nonlinear 
cubic array is given by 

A(3) = - i^cvc„ « " J uo + 2u§y cos(fc). (11) 

In comparison with the previously obtained result for 
the PN barrier ^21], Eq. (|11|1 adds an extra, albeit small, 
term that shows how the PN barrier is modified for the 
mode initially kicked in the lattice. Indeed, besides the 
first term dependent on the mode amplitude, Eq. (|ll|l in- 
cludes a linear term proportional to the factor y cos(fc), 
whose magnitude could be modified by a judicious ad- 
justment of the waveguide couplings and/or the value of 
the initial kick. 

Dependence of the PN barrier on the mode coupling 
suggests that, if we wish to find a way to engineer the 
value of the PN barrier in the lattice, we should study the 
properties of a modified model described by the evolution 
equation 

j^T^ + K+iUn+i + K-iw„-i + 7|u„pu„ = 0, (12) 
az 

where the coupling Vn between two neighboring guides 
is assumed to vary either through the effective propaga- 
tion constant or by a change in the spacing between the 
neighboring waveguides. To study the beam steering in 
this novel model, we use again as an initial condition the 
sech-like profile ^ , although this is not really fundamen- 
tal limitation, as argued below. 

We mention that a variation of the waveguide coupling 
in the array constitutes the starting point for our concept 
of the waveguide array engineering. A change of the cou- 
plings breaks the symmetry between the beam motion to 
the right and left at the moment of trapping, thus elimi- 
nating chaotic trapping observed in the case of homoge- 
neous arrays. 

We have tested different types of modulation in the 
array coupling and the corresponding structures of the 
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FIG. 3: Digitized controlled switching of a discrete soliton 
in a cubic nonlinear waveguide array by varying the beam 
intensity, for a fixed angle (fc — —0.9). Inset: Optimized 
modulation of the couplings Vn- 
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FIG. 4: Same as in Fig.|21 but in the engineered waveguide ar- 
ray with the coupling modulation shown in the inset of Fig. 13 

waveguide super-lattices. An example of one of such op- 
timized structure, where we modulate the coupling pa- 
rameter Vn in a step- like manner, is shown in the inset of 
Fig. O which also shows the discrete position of the soli- 
ton at the output as a function of the amplitude of the 
input beam, at a fixed value of the steering parameter 
k = —0.9. In a remarkable contrast with other studies 
(see, e.g., Ref. 0), the couphng modulation allows to 
achieve a controllable digitized switching of discrete opti- 
cal solitons in the array with very little or no distortion. 

As is shown in Fig. ^ by decreasing the amplitude 
of the input pulse at a fixed value of the steering pa- 
rameter k (in our example fixed to be fc = —0.9), it is 
possible to achieve self-trapping of the discrete soliton 
by the lattice at some (short) distance from the input at 
different waveguide positions. Due to the step-like mod- 
ulated coupling, we create a selection between the beam 
motion to the right and left at the moment of trapping 
thus suppressing or eliminating the chaotic trapping ob- 
served in homogeneous waveguide arrays. In this way, we 
achieve a controllable digitized nonlinear switching where 
the continuous change of the amplitude of the input beam 
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FIG. 5: Digitized controlled switching of a cubic discrete soli- 
ton by varying the beam input angle, for a fixed beam inten- 
sity. Inset: Optimized modulation of the couplings Vn- 



results in a quantized systematic displacement of the out- 
put beam by an integer number of waveguides. Conse- 
quently, for the parameters discussed above we observe 
almost undistorted switching up to eleven waveguides. 
Incidentally, we notice here that the use of a linear ramp 
potential (e.g., in the form Vn — an) for this purpose does 
not lead to an effective switching but, instead, it makes 
the soliton switching even more chaotic due to the phe- 
nomenon of Bloch oscillations which become randomized 
in the nonlinear regime. 

In Fig. 5 we show another example of the optimized 
coupling modulation, this time as a function of the ef- 
fective input 'kick', for a fixed beam intensity. In this 
case, we can achieve completely controlled switching up 
to nine waveguides. 

If the input beam was to excite initially five waveguides 
instead of three creating in this manner a wider excita- 
tion (i.e., being closer to the continuum limit), one could 
expect a smaller amount of radiation emitted. However, 
this would imply a longer distance before the beam gets 
trapped by one of the waveguides in the array due to 
the effective PN potential. Also, this means that one 
could, in principle, switch the soliton beam to any de- 
sired waveguide in the waveguide array, no matter how 
far; it would be just a matter of choosing an initial beam 
wide enough, i.e., closer to the continuum (in addition 
to optimize the coupling in a step- wise manner), by re- 
moving the random selection between the directions and 
suppressing the beam random switching. 

Another observation is that the scch-likc initial pro- 
file is not really fundamental. We have verified that the 
similar dynamics is observed for other types of the input 
beam profiles, including a nonlinear impurity-like input 
of the form 



w„(0) = A 



1 + 



2 \ |ri-nc|/2 



-ik{n—nc) 



(13) 



the observation that for any system with local nonlin- 
earity a narrow initial profile will render the system into 
an effective linear one containing a small nonlinear clus- 
ter (or even a single site); the bound state will therefore 
strongly resemble that corresponding to a nonlinear im- 
purity |2^. 



B. Continuous model 

In Sec. Ill A1 we have obtained exceedingly interesting 
results for the switching of discrete solitons, via the use 
of the discrete model and tight-binding approximation. 
In this section, in order to confirm our predictions, we 
perform the corresponding numerical simulations of the 
continuous evolution equations of the electric field in- 
side a waveguide array, by means of the Beam Propa- 
gation Method (BPM) '23|. As we show below, these 
results support our major findings concerning the digi- 
tized switching of an optical beam in engineered nonlin- 
ear waveguide arrays. 

In the continuous model, the starting point is an op- 
tical beam propagating in a three-dimensional medium 
with the refractive index that varies in space and is also 
intensity-dependent accounting for the Kerr effect. For 
a nonlinear waveguide array, the beam is assumed to 
propagate along the z direction and to diffract or self- 
focus along the transversal directions x and y. Assum- 
ing that the beam envelope A(x^ y, z) varies with z on a 
scale much longer that the wavelength A, the beam enve- 
lope is found to obey [23| the three-dimensional nonlinear 
Schrodingcr (NLS) equation. 



r. n dA 
oz 



d'^A &^A 



dx^ dy^ 



2/3ofcon2 
S 



\A\'^A^O (14) 



where /3o = 27rno/A is the beam propagation constant, 
n2 is the Kerr coefficient of the nonlinear guide, S is the 
area of the mode, and \A\'^ is the beam power. 

The waveguide structure creates a periodic modula- 
tion of the refractive index in only one of the transversal 
directions (say, x). Along the other transversal direc- 
tion y we assume the beam confinement. Therefore, the 
real dimensionality of the system is two (i.e., the lon- 
gitudinal propagation and the transversal spreading or 
self-trapping) . We write the electric field envelope in the 
form A( x, y , z) = a(x, z)b{y) and use the effective index 
method |23| to formally transform Eq. (|14l) into an effec- 
tive two-dimensional equation 



1 d^a fcon2 



.da 



dz 2kon^Q{x) dx^ A 



\a\^a = 0, (15) 



The reason for this universal behavior seems to rest on 



eff 

where A^q — S \b\'^dy/ \b\'^dy, is the effective 
area over which the nonlinear interaction occurs. Pa- 
rameter rig-gp is the effective, space-varying linear index 
of refraction for the one-dimensional problem. For our 
problem, nggr(a;) consists of a periodic array of parallel 
slabs with the indices uq and ng -I- An. 
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FIG. 6: Results of numerical BPM simulations of the continuous model for the switching of a 1.5/^m Gaussian beam propagating 
in an array of 41 AlGaAs waveguides as a function of the input beam intensity, for a fixed input beam angle (0.86°). Left: 
2200 W, center: 2330 W, right: 2700 W. 



We solve Eq. H15|) numerically by the well-known beam 
propagation method (BPM), using typical experimen- 
tal parameter values. In our simulations, we use an 
array of 41 waveguides, 10mm long, a Kerr coefficient 
712 = 1-5 X 10~^^TO^/iy, assuming a Gaussian beam of 
the height 3/im, 8/im wide, at A = 1.5/^77i, and 0.86° 
input angle. The effective one-dimensional medium con- 
sists of a periodic array of slabs, 4/ir7i wide, with 6fim 
center-to-center separation, uq = 3.27, with a modula- 
tion of the refractive index An = 0.0014. 

Figure shows the beam switching as the initial input 
beam power is varied: Switching to seven, five and three 
guides was observed for a power of 2200 W, 2330 W and 
2700 W, respectively. As expected from the estimates 
of section III A 21 as power is decreased, so does the PN 
barrier, which allows the beam to get self-trapped farther 
away from the vicinity of the input guide. 

The above constitute a strong confirmation of the ro- 
bustness of the switching mechanism predicted by a sim- 
ple theory, against more realistic effects. Fine tuning of 
the switching by an appropriate modulation of the effec- 
tive waveguide couplings will allow a fine control of the 
position of the destination waveguide. However, this task 
is beyond the main scope of the present paper. 



phase modulation of the fundamental as well as the sec- 
ond harmonic waves p^ . This phase modulation ac- 
companies the familiar amplitude modulation, being the 
basis of any frequency conversion, and it may produce 
the effects which resemble those known to occur in cu- 
bic nonlinear materials. Typical examples are all-optical 
switching phenomena in intcrfcrometric or coupler con- 
figurations as well as the formation of spatial and tempo- 
ral solitons in planar waveguides (see, e.g., Ref. and 
the references therein). 

Recently, it was demonstrated theoretically p^l23.l2^ 
that arrays of quadratic nonlinear waveguides represent a 
convenient system to verify experimentally many theoret- 
ical predictions for the dynamics of nonlinear latices with 
cubic nonlinearity. The first experimental observation of 
discrete quadratic solitons has been reported recently by 
Stegeman and co-authors , who demonstrated the for- 
mation of discrete quadratic solitons in periodically poled 
Lithium Niobate waveguide arrays, excited with funda- 
mental wave pulses at a wavelength of 1572 nm. These 
experimental observations open many perspectives for 
employing much larger nonlinearities provided by non- 
linear quadratic materials. In this section, we extend the 
concept of the controlled digitized soliton switching dis- 
cussed above to the case of quadratic discrete solitons. 



III. QUADRATIC NONLINEAR WAVEGUIDES 

Up to now we discussed the arrays of weakly coupled 
waveguides with the cubic nonlinearity. However, during 
last years a growing interest is observed in the study of 
nonlinear optical effects based on the so-called quadratic 
nonlinearities. In contrast to the conventional studies 
of quadratic nonlinearities where the main attention is 
centered primarily on parametric processes and the fre- 
quency conversion, more recent works are focused on the 



A. Discrete Model 

The standard discrete model for an array of weakly 
coupled quadratic nonlinear waveguides has the form |3] : 

+ K(an+i + an-i) + 2726„a* = 
i^ + Vbibn+i+bn-i)+(3bn+l2al = 0, (16) 
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where a„ and 6„ represent the amphtudes for the funda- 
mental ((jj) and second harmonic {2lo) fields in the n-th 
guide, Va and VJ, stand for the linear couplings between 
the nearest-neighbor waveguides. Parameter 72 describes 
the nonlinear second-order coefficient proportional to the 
second-order dielectric susceptibility, and (3 is the effec- 
tive mismatch between the fields in the array. 

As in the case of the cubic nonlinearity, the system 
possesses two conserved quantities: the total power, 



P 



E 



(|a„(z)|2 + 2|6„(z)|^ 



(17) 



and the system Hamiltonian, 



H = -^(K<a„+i+F,fe;6„+i + (/3/2)|6„ 

n 

+ l2alb*^ + c.c.). 



(18) 



However, unlike the case of the cubic nonlinear waveg- 
uide arrays where it is possible to find analytical solutions 
in the continuum limit which can be used as input pro- 
files for numerical simulations of discrete solitons, in the 
case of the quadratic nonlinearities no exact solutions are 
available. Thus, we should resort to the limit of strongly 
localized modes (SLMs) in order to calculate the PN bar- 
rier and use the SLM profile as an input beam profile for 
the numerical computation of the soliton switching. 



B. Localized modes and the PN barrier 

As in the case of the cubic nonlinearity, we calculate 
the PN barrier as a difference between the values of the 
Hamiltonian for the odd and even strongly localized two- 
component modes. 

Odd modes. We search for approximate solutions of 
Eq. Hlt)|) of the form 

an = ao{....,0,ale*^l,ale-'^0,...} e^^i^ 

bn - 6o{....,0,6le2*^l,6le-2*^0,...} e2^^l^ (19) 

where Oq and 60 the amplitudes of two harmonics 
composing a localized mode, k is the initial beam angle 
or effective 'kick', Ai is the longitudinal propagation con- 
stant, and ai and 61 are small parameters that should be 
determined from the equations. After substituting the 
ansatz H19|l into Eqs. (|16|l and keeping only linear terms 
in oi and hi, we obtain: Ai — 2aiVa cos(fc) + 27269, Uq = 
(&o/72)[2Ai-/3-2feiVi,cos(2fc)], ai = (Va/Ai) cos(fc), and 
bi — Vb cos(2fc)/(2Ai — (3). From these relations, we find 
Ai w 272&0J which implies 



46^ 



ai 



- {f3/l2)b0: 

Vg cos(fc) 
27260 



< 1, bi 



Vb cos(2fc) 
47260 - (3 



< 1. (20) 



Even modes. Now we search for approximate solutions 
of Eqs. (|16|l of the form 



^„ = ao{••,0,ale'^l,e '''^aie 2'^0, ...}e 



1X2 z 



6o{--,0,6ie 



2ik 



l,e 



-2ik 



6le"**^0,...}e2'^^^ (21) 



where fio and 60 are the amplitudes of the coupled har- 
monics, k is the initial beam angle or effective 'kick', 
A2 is the longitudinal propagation constant, oi and bi 
are small parameters determined from the equations 
of motion. After substituting Eq. into Eq. l(TB|) 

and keeping only linear terms in di and 61, we obtain: 
A2 = (1 + ai)Kcos(fc) + 27260, al = (6o/A2)(2A2 - 
/3 - (1 + 6i)Hcos(2fc)), 5,1 = (K/A2)cos(fc) and 61 = 
VbCOs(2fc)/(2A2 — /?). From these relations, we find 
A2 ~ Va cos{k) + 27260, which implies 



^.2 



0.1 



61 



46g - (60/72) [/? - 214 cos(fc) + Vb cos{2k)], 

Va COs(fc) 



Va COs(fc) + 27260 

Vb cos(2fc) 



« 1, 



< 1. 



(22) 



2Va cos(fc) + 47260 - (3 

From Eq. I|20() and Eq. 122() . it is easy to see that, in or- 
der to obtain SLM, the nonlinear term 7260(60) should be 
be much larger than the linear coupling terms, Va cos{k) 
and Vb cos(2A;). 

From Eq. (|17|1 we calculate the total power 



Podd»ao + 26^ + 0(af,6f), 



(23) 



Pevc„«2a2 + 462 + 0(^2, 6?), (24) 

and the Hamiltonian of each mode, 

Hodd ~ —SalaiVa cos{k) — Sb^biVb cos{2k) 

~pbl-2j2albo + 0{al,bl), (25) 

i^cvcn ~ -45^(1 + 2hi)Va COs(fc) - 472^^60 

-46g(l + 2bi)Vb cos(2A:) - 2/36^ + 0(&?, 6?). (26) 

We follow the same reasoning as in the case of the cu- 
bic nonlinear waveguide arrays and calculate the effective 
PN barrier. Such calculations look simpler for the physi- 
cally important case of vanishing mismatch, /3 « 0. With 
that assumption, and imposing that the power content of 
both, odd and even modes, are equal, Podd = -Pcvcn, we 
obtain 36g w 66g + (60/72) ( 2Va cos(fc) - Vb cos(2fc) ), and 
then 

~ 60 ( 2Va cos(/c) - Vb cos(2/c) ) 

'''"71 ^2 • ^''^ 

In terms of 60, the Hamiltonian of both the modes can 
be approximated as 



Hodd - -Sj2bt + Oiai,bi) 



(28) 
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FIG. 8: Switching to 6, 5 and 4 sites of a discrete quadratic 
soliton (SLM) with a slight intensity variation of the funda- 
mental mode. 



FIG. 7: Controlled digitized switching of a discrete quadratic 
soliton by a variation of the intensity of the fundamental 
mode, for a fixed input 'kick'. Inset; Optimized modulation 

of Vn = Va,n = Vb,„. 

TJeven ~ -4^2725}^ " 8^6^ COs(fc) 

-2Vbbl cos(2fc) + 0{ai • ai, 61 • 61). (29) 

Finally, we calculate, in this approximation, the PN bar- 
rier of the strongly localized modes, 

A(2) - H a,- H w 

-8C72fcg + 2hl{ Wa cos(fc) -f Vb cos(2fc) ), (30) 

where C = (1 - \/2/2). The PN barrier of an array of 
nonlinear quadratic waveguides (|30|l has been obtained, 
to the best of our knowledge, for the first time. It shows 
some interesting features: The main term l|5n)l is cubic in 
the mode amplitude, while for the cubic case it was quar- 
tic [see Eq. Also we notice that the first correction 

to the PN barrier (|3(J|I is linear in the couplings, and it 
depends on the square of the SLM amplitude. This is ex- 
actly the same term as in the case of the nonlinear cubic 
array. This implies that the first-order correction is more 
important in the nonlinear quadratic array than that in 
the nonlinear cubic array suggesting that the appropri- 
ate engineering of the couplings and/or input 'kick' to 
achieve digitized switching should be easier to achieve. 

For the numerical simulations, we use the initial input 
in the form of an odd mode, 

a„(0) = ao•al"-"^le-'("-"=)^ 

h,m = 6o•6^"-"^le-2'("-"=)^ (31) 

for n — Tie = 0, ±1, and a,i(0) — 5„(0) = 0, otherwise. 
In Eq. (EIJ we use ao ~ \/46g - (/3/72)&o, ai ~ K/272&0 
and 61 « 14/(472^0 ~ /3)- We consider an array of 41 
waveguides with the initial input centered at the mid- 
dle, ric = 20. For simplicity, we also assume the case 
of complete phase-matching, i.e. /3 « 0, and identical 
coupling for both the harmonic fields, Va,n = 14, «• Fig- 
ure [7| shows that the controlled digitized switching of the 



discrete two-frequency (fundamental -I- second harmonic) 
soliton can be achieved for up to 10 waveguides, by vary- 
ing the intensity of the input fundamental mode for a 
fixed parameter fc = — 1. The inset of Fig. [7| shows the 
coupling modulation required to achieve this type of en- 
gineered soliton switching, which is particularly simple 
and consists of only a single change of about 5% in the 
value of the coupling parameter. Figure |H1 demonstrates 
the switching of the discrete mode, composed of the fun- 
damental and second- harmonic fields, to six, five, and 
four neighboring waveguides, as the intensity of the in- 
put fundamental mode is increased. In this respect, it 
is interesting to point out that in all cases of the digi- 
tal switching both the fundamental and second-harmonic 
fields act as a strongly coupled state, and that no 'lagging 
behind' was observed of any of the modes with respect 
to the other. 

We have performed other simulations with the 
quadratic nonlinear array including the cases 14 = (de- 
coupled second-harmonic fields in the array) and VJ, — 
aVa (reduced coupling of the second-harmonic fields) 
with a < 1, etc. In all of those cases, we have observed 
the digitized switching of the discrete solitons by engi- 
neering the coupling in the array as discussed above. 



IV. CONCLUSIONS 

We have suggested and demonstrated numerically a 
simple but yet effective method for controlling nonlinear 
switching of discrete solitons in arrays of weakly cou- 
pled optical waveguides. We have demonstrated how to 
achieve the digitized switching of discrete optical solitons 
in weakly coupled arrays of cubic and quadratic nonlin- 
ear waveguides described, in the framework of the tight- 
binding approximation, by discrete models such as the 
DNLS equation with a step-like variation of the waveg- 
uide coupling parameter. Our approach involves a weak 
step-like modulation of the coupling strength (or, equiv- 
alently, distance between the waveguides) in the arrays 
with the period larger than the waveguide spacing. Such 
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kind of a super-lattice waveguide structure allows to mod- 
ify the trapping properties of the array due to discrete- 
ness as well as engineer the strength of the effective 
Peierls-Nabarro potential arising due to the lattice dis- 
creteness. In particular, we have demonstrated the digi- 
tized switching of a narrow input beam for up to eleven 
waveguides, in the case of the cubic nonlinear array, and 
up to ten waveguides, in the case of quadratic nonlinear 
array. We have confirmed our predictions for a full-scaled 
continuous model and realistic parameters by employing 
the beam propagation method. 
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